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ABSTRACT 

Spectroscopic redshifts have been measured for 9 /f-band luminous galaxies at 1.7 < z < 2.3, selected with 
Ks < 20 in the K20 survey region of the Great Observatories Origins Deep Survey area. Star formation rates 
(SFRs) of ~ 100-500 Mq yr are derived when dust extinction is taken into account. The fitting of their 
multi-color spectral energy distributions indicates stellar masses M £ 10 M© for most of the galaxies. Their 
rest-frame UV morphology is highly irregular, suggesting that merging-driven starbursts are going on in these 
galaxies. Morphologies tend to be more compact in the near-IR, a hint for the possible presence of older stellar 
populations. Such galaxies are strongly clustered, with 7 out of 9 belonging to redshift spikes, which indicates 
a correlation length ro ~ 9-17 h~ l Mpc (la range). Current semianalytical models of galaxy formation appear 
to underpredict by a large factor ( ^ 30) the number density of such a population of massive and powerful 
starburst galaxies at z ~ 2. The high masses and SFRs together with the strong clustering suggest that at z ~ 2 
we may have started to explore the major formation epoch of massive early-type galaxies. 
Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift — galaxies: starbursts 
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1. INTRODUCTION 

The remarkable success of the Cold Dark Matter (CDM) 
scenario to account for the cosmic microwave background 
power spectrum (Bennett et al. 2003), leaves understand- 
ing galaxy formation and evolution as one of the most com- 
pelling, unresolved issues of present cosmology. Semiana- 
lytical renditions of CDM hierarchical paradigm have so far 
favored a slow growth with time, with a major fraction of the 
mass assembly taking place at z 1 (e.g., Baugh et al. 2002; 
Somerville, Primack, & Faber 2001), with virtually all mas- 
sive galaxies disappearing by z ~ 1.5. Recent results from 
the K20 project (Cimatti et al. 2002a,b,c; Daddi et al. 2002; 
Pozzetti et al. 2003) appear to be at variance with these expec- 
tations. The K20 project consists of a spectroscopic survey of 
~ 500 Ks < 20 objects selected over 52 arcmin 2 , and has re- 
vealed a sizable high redshift tail in the galaxy redshift distri- 
bution, where ~ 30 galaxies (~ 6% of the total sample) were 
found at z > 1.7. Semianalytical CDM models would have 
predicted no galaxy at all at such high redshifts in the whole 
sample (see Fig. 4 in Cimatti et al. 2002c). The redshift 
distribution could be reproduced with pure luminosity evolu- 
tion (PLE) models, although not for all realizations (see also 
Somerville et al. 2003). However, for only a few among the 
Z ^, 1.5 galaxies was a spectroscopic redshift available, while 
for all other such galaxies only the photometric redshifts 
could be obtained. In order to put on firmer grounds such 
major result of the K20 project (and to understand the nature 
of these high-z galaxies) we have conducted new VLT spec- 
troscopic observations of galaxies with either photometric or 

1 Based on observations collected at the European Southern Observatory, 
Chile (ESO programs 70.A-0140, 168.A-0485), and with the NASA/ESA 
Hubble Space Telescope, obtained at the Space Telescope Science Institute, 
which is operated by AURA Inc, under NASA contract NAS 5-26555. 

2 ESO, Karl-Schwarzschild-Str. 2, D-85748 Garching, Germany 

3 Osservatorio Astrofisico di Arcetri, L.go E. Fermi 5, Firenze, Italy 

4 Caltech, Mail code 105-24, Pasadena (CA) 91 125 

5 Osservatorio Astronomico di Bologna, via Ranzani 1, Bologna, Italy 

6 Osservatorio Astronomico di Trieste, via Tiepolo 1 1 , Trieste, Italy 

7 Racah Institute for Physics, The Hebrew University, Jerusalem, Israel 

8 Oss. Astron. di Roma, via DeH'Osservatorio 2, Monteporzio, Italy 



(uncertain) spectroscopic redshift above z ~ 1.7. This let- 
ter reports the results of the new spectroscopic observations, 
and combines them with the optical HST+ACS and infrared 
VLT+ISAAC imaging made available by the Great Observa- 
tories Origins Deep Survey (GOODS) project (Giavalisco et 
al. 2003). We assume a Salpeter IMF; A , n M = 0.7,0.3, and 
h = // [km s" 1 Mpc^l/lOO = 0.7. 

2. THE SPECTROSCOPIC SAMPLE 

A sample of 20 galaxies with photometric redshifts 
Zphot ~ 1 .7 were selected among the 41 galaxies without spec- 
troscopic redshift identification in the 32 arcmin 2 K20 field 
that is included in the GOODS-South area. The photometric 
redshifts where improved over an earlier estimate (Cimatti et 
al. 2002b) by including the ultra-deep JHKs photometry from 
the GOODS VLT+ISAAC imaging. Within this sample, 10 
objects with most secure z p hot have been observed in Novem- 
ber 2002 with VLT+FORS2, integrating for 7.2ks with 0".6 
seeing, and using the 300V grism covering the range 3600- 
8000 A with 13 A resolution for a 1" slit. 

The spectra were reduced and calibrated in a standard way 
(Cimatti et al. 2003b) and co-added to already existing spec- 
tra, when available. Redshifts were measured for 7 galax- 
ies from absorption features in their blue continua identified 
as UV metal lines (Fig. [Q. For the z > 2 galaxies Lya 
in absorption and the onset of Lya forest are also detected. 
One galaxy (ID#5) shows HeIIA1640 and CIII]A1909 emis- 
sion lines. Hints for those emission lines are found also for 
object ID#9, for which redshift is less secure because of the 
faint and noisy spectrum. Together with 2 previously iden- 
tified galaxies (ID#1 and ID#2), a sample of 9 galaxies with 
spectroscopic redshift z spec > 1 .7 is now available among the 
304 Ks < 20 galaxies in the K20/GOODS-South area (Ta- 
ble 1). Correspondingly, the fraction of Ks < 20 galaxies 
with z spec > 1.7 is 3.0! 2 -^%, for a surface density of f f 



arcmin 



and a comoving density of 4.6+^ x 10" 



h 3 Mpc" 

(the range 1.7 < z < 2.25 is used, hereinafter, for volume 
calculations). The uncertainties are poissonian at the 95% 
c.l. Accounting for cosmic variance due to clustering (Sect. 
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FIG. 1. — The panels show the spectra of, from top to bottom: (1) 
ID#7 at z = 2.227; (2) ID#5, the X-ray and radio source at £ = 2.223; 
(3) the average spectrum of 8 z > 1.7 K-band luminous galaxies (ID#2 
has only a near-IR spectrum showing Ha emission) (4) the average spec- 
trum of ~ 250 LBGs with Lya in absorption (Shapley et al. 2003). The 
principal metal lines observed in the UV for starburst galaxies are marked 
for reference, namely the ISM lines SiIIA1260; 1304; 1527; 1534;, OIA1303, 
CIIA1334.5, CTVA1548; 1551, SiIVA1394; 1403, FeIIA1608, A1IIA1671, 
A1IIIA1855; 1863; as well as the photospheric absorption lines (shown with 
an asterisk) FeIIIA1926, SVA1502, SiIIIA1294, SiIIA1485,CIIIA1427. 



[3ji would significantly alter only the upper bounds (see e.g. 
Eq. 8 in Daddi et al. 2000). These densities would increase 
by a factor of two by including the remaining objects with 
Zphot £1.7. The good agreement between z pn ot and z spec in 
the poorly tested region 1.7 < z < 2.3 (Table 1) suggests a 
fair fraction of the latter to be genuine z ~ 2 galaxies. Ob- 
taining z S pec for these is difficult due to their redder colors and 
thus fainter optical magnitudes with respect to galaxies with 
measured z spe c. 

3. PROPERTIES OF K-BAND LUMINOUS GALAXIES AT z ~ 2 

Using the wealth of subsidiary data available on the 
GOODS-South area we investigate in this section the nature 
of these /T-band luminous galaxies with z spe c ~ 2. 
- Star formation rates. The rest- frame UV spectra of these 
galaxies indicate that they are actively star-forming. They are 
very similar to the template Lyman Break Galaxies (LBG) 
spectrum (Shapley et al. 2003), and often show SUIIA1294 
absorption due to OB stars (Fig. Lya emission is never 
detected, as for 25-50% of LBGs (Shapley et al 2003), a hint 
for significant dust extinction. The UV fluxes at 2800A corre- 
spond to star formation rates (SFRs) in the range of ~ 10^-0 
M yr before extinction correction (Kennicutt 1998). The 
extinction at 1600 A was estimated from the UV spectral slope 
(Meurer et al. 1999) derived from the multicolor photom- 
etry. The SFRs derived in this way (adopting the Calzetti 



et al. 2000 extinction law) are very high, with a median of 
« 400 Moyr -1 . As an alternative estimate hyperz was used 
(Bolzonella et al. 2000) to fit Bruzual & Chariot model spec- 
tral energy distributions (SED) to the observed ones from the 
available deep VLT BVRIzJHKs photometry, assuming con- 
stant SFR (CSF), and solar metallicity. The best fitting models 
require reddening in the range E(B-V) ~ 0.3-0.6 and intense 
SFRs up to ~ 500 M yr _1 , with a median of 150 M yr _1 . 
Using a SMC extinction law would lower both the SFRs and 
E(B-V) estimates. Interestingly, ID#5 (which has the high- 
est extinction-corrected SFR) is a faint (soft) X-ray source 
(XID563 in the Chandra Deep Field South catalog, Giacconi 
et al. 2002). If completely due to star formation, its X-ray 
luminosity i^-iokev ~ ^-7 X 10 42 erg s -1 cm -2 corresponds to 
SFR « 500 M yr _1 (Nandra et al. 2002). The object is also a 
faint 1.4 Ghz radio source with a flux density of 103 ± 13 /iJy 
(Kellermann et al. 2003, private comunication), fully consis- 
tent with the tight X-ray-radio luminosities correlation shown 
by Ranalli et al. (2003) for actively star-forming galaxies. We 
cannot definitively rule out the presence of an obscured AGN, 
but the lack of AGN signatures in its UV spectrum, showing 
instead a strong SHIIA1294 photospheric absorption line, the 
non detection in the Chandra hard band and the low X-ray 
to optical flux ratio (/og(/o.5-2kev/./}0 ~ -1-3) indicate this is 
a vigorous star-forming galaxy. The stacked X-ray fluxes of 
the undetected sources give a 2a detection corresponding to 
< SFR > pa 100 M yr (r = 2.1 is assumed following Brusa 
et al. 2002). In general the limits on the X-ray luminosities 
(Lx 10 42 erg s -1 cm -2 ) and the low X-ray to optical flux 
ratios (log(fo.s-2kev / /r) ~ - 1-5) imply that our sample con- 
tains virtually no AGN. The high SFRs qualify these galaxies 
as starbursts, and allow to build up the equivalent of a local 
M* ~ lO n M galaxy in 0.2-1 Gyr. 

- Stellar masses. Although the SEDs are reddened in the 
rest-frame UV, they appear even redder toward the near-IR, 
where they show a steep flux increase, starting in the F850LP 
band or beyond, which is suggestive of relatively old stellar 
populations. The CSF models discussed above imply M stars = 
0.3-5.5 x 1O U M , and luminosity-weighted ages of about 
250-1700 Myr. In this case, the near-IR bump is reproduced 
by a prominent Balmer break. The CSF models are likely 
to underestimate the masses of the galaxies, as older stars, 
with higher mass to light ratios, may well be present and yet 
their light would be outshined by the younger ones. In order 
to estimate reliable upper limits, the minimal contribution to 
the A"-band light by the ongoing starburst is determined by 
fitting a very young ( 10 Myr) reddened component to the 
SEDs between the B and / bands. This component accounts 
for 30-50% of the /T-band light. Assuming the remaining K- 
band light is due to a maximally old 3 Gyr stellar-population 
component, the resulting masses are typically a factor of 2-5 
higher than estimated from CSF models, similarly to what is 
found for LBGs (Papovich et al. 200 1 ) . 

- Clustering. Significant redshift pairing is observed 
among z ~ 2 galaxies, a clear indication of strong cluster- 
ing. Monte Carlo simulations are used to constrain the corre- 
lation length ro from the short scale pairing, assuming a slope 
7 = 1.8 for the correlation function (Daddi et al. 2002). A flat 
selection function between z = 1 .7 and z = 2.25 is used. Seven 
independent pairs within 5 h~ l comoving Mpc are found in 
our sample of 9 galaxies, implying ro > 7 h Mpc comoving 
(95% c.l.), and a most likely range 9-17 /T 1 Mpc (68% c.l.). 



Near-IR bright galaxies at z — 2 
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FIG. 2. — ACS (top and center row for F435W and F850LP respectively, epochs 1+2+3) and VLT+ISAAC (bottom row for Ks, seeing 0.5 ) imaging for the 
galaxies with spectroscopic identification. The images are 5" on a side. Redshift measurement for each galaxy is given in the bottom panels. 



Table 1. a:-band luminous starbursts at z > 1.7 in the K20/GOODS area 
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1 


1.727 


1.74 


19.99 


1.44 


3.15 


-0.7 


26 


93 


0.3 


0.36 


3.4 


0.6 


4.8 


2.2 


0.21 


0.9 


2 


1.729 


2.54 


19.07 


1.99 


4.54 


0.7 


13 


490 


0.6 


0.25 


11 


0.8 


7.1 


1.8 


0.30 


1.4 


3 


1.901 


1.65 


19.68 


1.57 


3.42 


-0.9 


27 


155 


0.3 


0.52 


7.9 


0.6 


5.4 


1.9 


0.41 


0.6 


4 


2.060 


1.78 


19.31 


1.83 


4.15 


-0.2 


35 


487 


0.5 


0.50 


25 


0.7 


5.9 


2.0 


0.34 


1.3 


5 


2.223 


2.40 


18.72 


2.15 


4.15 


0.0 


38 


540 


0.4 


1.0 


55 


0.4 


3.5 


3.2 


0.25 


0.3 


6 


2.226 


2.24 


19.94 


1.69 


3.30 


-0.4 


18 


88 


0.3 


0.71 


6.0 


0.7 


5.5 


2.6 


0.33 


0.5 


7 


2.227 


2.11 


19.45 


1.80 


3.46 


-1.2 


31 


178 


0.3 


0.72 


13 


0.7 


5.5 


2.2 


0.41 


1.3 


8 


2.228 


2.43 


19.74 


2.02 


3.66 


-0.4 


26 


121 


0.3 


1.4 


17 


0.9 


7.8 


2.1 


0.25 


0.9 


9 


2.25 


2.29 


19.94 


1.94 


3.65 


-0.7 


30 


153 


0.3 


1.7 


26 


1.1 


9.4 


2.7 


0.29 


1.2 



NOTE. — (1) The redshift for ID#9 is less secure. (2) UV spectral slope. Typical errors are ±0.1. (3) SFR derived from the 2800A luminosity without extinction correction. (4) 
Star formation rate, extinction, luminosity weighted stellar age and stellar mass derived from SED fitting of CSF models with reddening. (5) Quantities measured in the F850LP band, 
typical errors are AC, AA, AS = 0. 15,0.15,0.05. 



- Morphology. In the HST+ACS and VLT+ISAAC images 
taken for the GOODS project all galaxies show a rather ir- 
regular light distribution (Fig. |2j, with bright knots and low 
surface brightness regions, often split into separated compo- 
nents. We measured the CAS parameters (Conselice 2003; 
Bershady et al. 2000; and references therein), finding rela- 
tively high dumpiness (S), high asymmetry (A) and very low 
concentration (C) (see Table [Q. As S is known to correlate 
with the SFRs, the large S values are consistent with the high 
SFRs estimated above. The A values of most galaxies are 
consistent within the errors with the limit of A > 0.35, typical 
of galaxies undergoing merging or that experienced merging 
in the last Gyr (Conselice 2003). The low C values are also 
typical of local merging-driven starbursts, or ultra luminous 
infrared galaxies (ULIRGs). There is a trend for increasing 
C from the rest frame far-UV (F435W band) to the optical 
(/£-band, resolution effects having been taken into account), 
implying a morphological K-correction. Also this is typical 
of starburst galaxies (see e.g. Dey et al. 1999; Smail et al. 
2003), and may indicate the presence of an older bulge/disk 
component (e.g. Labbe et al. 2003b) or a higher reddening in 
the central regions. All the galaxies appear rather extended, 
allowing to host the high estimated SFRs. The average half 
light radius is ra — 0".7 in the F850LP band, about ~ 6 kpc. 

4. RELATING TO OTHER z > 2 GALAXY POPULATIONS 

We now compare the properties of these /T-band lumi- 
nous galaxies to those of other relevant populations at z ^ 2, 
namely: LBGs at z ~ 3, very red z > 2 galaxies, and SCUBA 



sources. Compared to LBGs (e.g. Giavalisco et al. 2002), 
these near-IR bright starbursts at z ~ 2 have, on average, larger 
sizes, higher masses and SFRs, and stronger clustering. De- 
spite their spectral similarity, these galaxies are not just a spe- 
cial subsample of LBGs. Indeed, it appears that they have 
redder UV continuum than the reddest LBGs template (Fig. 
[l). In fact, most objects in Table 1 have (3 > -1, while most 
LBGs have UV slopes between (3 = -2 to -1, and virtually 
none has (3 > -0.5 (Adelberger & Steidel 2000). Hence, the 
two populations appear only partially overlapping. 

These z ~ 2 starbursts are red in the near-IR, with J — 
K 1.7, and their clustering is consistent with that of much 
fainter Ks < 24 galaxies at z > 2 with J-K > 1.7 colors 
(Daddi et al. 2003). In fact, van Dokkum et al. (2003) 
found significant redshift pairing among 5 galaxies at z > 2 
selected with J-K > 2.3 (Franx et al. 2003; Labbe et al. 
2003a). Nevertheless, the two samples show different proper- 
ties, as strong Lya emissions, regular morphologies, and AGN 
signatures are common among van Dokkum et al. objects. 
Our 9 spectroscopically confirmed galaxies have J-K < 2.3 
and very clumpy and asymmetric morphologies. We conclude 
that there is a large variety of properties among K-bmd bright 
galaxies at z > 2, that we are just starting to explore. 

Given the estimated SFRs, redshift range and peculiar mor- 
phology, some of our galaxies are potential SCUBA sources 
(Chapman 2003a,b). Red UV SEDs with (3 > -0.5 are in- 
deed common among SCUBA sources (e.g. Dey et al. 1999, 
Chapman et al. 2002), which also appear to have large clus- 
tering (e.g. Webb et al. 2003). Nevertheless, SCUBA sources 
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have much lower spatial density and are often much fainter 
than Ks = 20, while some of our sources may have too low 
SFRs to be submm-bright. Hence, also in this case the two 
populations are likely to overlap only partially. 

5. DISCUSSION 

These ,^-band luminous starbursts provide a substantial 
contribution to the cosmic SFR density (SFRD) at z ~ 2: 
adding up the SFRs from SED modeling we derive SFRD ~ 
0.04 M yr -1 Mpc -3 from the 9 spectroscopically confirmed 
galaxies alone. This estimate is certainly affected by in- 
completeness, and yet it already represents ~ 30-60% of the 
SFRD within the range 1.5 < z < 3 (see e.g. the compila- 
tion by Nandra et al. 2002). These galaxies are also among 
the most massive systems detected at z ~ 2. Six objects in our 
sample have M stars > 10 M© (conservative estimates), result- 
ing in a number density ~ 10~ 4 Mpc -3 and a mass density of 
~2x 1O 7 M Mpc" 3 , both w 10 % of the corresponding local 
value (Cole et al. 2001). Integrating over the mass function 
predicted by the Baugh et al. (2002) model at z = 1.92, one 
expects on average only 0.2 galaxies with M stars > 10 11 M Q 
within the explored volume, and even less if the semianalyt- 
ical mass function was properly normalized at z = 0. There- 
fore, these semianalytical models underestimate the number 
of massive galaxies at z ~ 2 by about a factor of 30, and pos- 
sibly much more given the incompletenesses of our spectro- 
scopic sample. The assembly of massive galaxies apparently 
took place at a significantly larger redshift (earlier epoch) than 
predicted by the models (see also Genzel et al. 2003). On the 
other hand, these z ~ 2 galaxies are too actively star-forming 
and irregular to be consistent with PLE models with high red- 
shift of formation. The agreement between the observed red- 
shift distribution at z > 1.7 and the PLE model described in 
Cimatti et al. (2002c) is therefore likely to be just chance. 

These /if -band luminous starbursts are very strongly 
clustered, suggesting they are hosted in very massive and 
biased environments, which itself argues for these objects 



being quite massive. At z < 2 the only known sources with 
r$^l h Mpc are old, passively evolving EROs (Daddi et 
al. 2000, 2001) and local massive ellipticals (Norberg et al. 
2002). These z ~ 2 galaxies are therefore likely to evolve 
into such classes of objects. If star-formation ends rapidly, 
it would take them ^ 1 Gyr to develop very red optical 
to near-IR colors and to morphologically relax to regular 
bulge-dominated galaxies. This scenario, with massive 
spheroids still forming at z ~ 2-3, would be quite in good 
agreement with some properties of z ~ 1 old EROs, including 
their number counts (Daddi, Cimatti & Renzini 2000), hints 
for residual star-formation present in their UV rest-frame 
(McCarthy et al. 2001), and with the inferred formation 
redshifts (2.4 ±0.3) of their stellar populations (Cimatti et 
al. 2002a). At the same time, it would predict a paucity of 
passive EROs at, say, z > 1.3-1.5. Finally, we notice that a 
major shift seems to happen for the clustering properties of 
star-forming galaxies from z ~ 1, where they have very low 
clustering (see e.g. Daddi et al. 2002), to z ~ 2, where they 
have a very large one. The straightforward interpretation is 
that while at z 1 star formation is mostly confined to low- 
mass galaxies, at z ~ 2 we are starting to see the major build 
up phase of massive early-type galaxies. It remains to be 
determined whether this z ~ 2 activity represents the peak or 
the low-z tail of the massive spheroid formation epoch. With 
the existing technology we should soon be able to answer this 
question, mapping massive galaxy assembly as a function of 
both redshift and large scale structure environment. 
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